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1. Introduction 



The definition of the double scaling limit was a significant breakthrough in physics 
because it provided a framework in which nonperturbative aspects of string theory could 
be discussed in a rigorous and precise way []1|-|3[] . Although the nonperturbative formulation 
given by the matrix models is not necessarily the "correct" one, it is the only one available 
and hence deserves close scrutiny. 

Unfortunately, the original hope that a Painleve I transcendent would provide the 
nonperturbative definition of the free energy of pure 2D quantum gravity evaporated with 
the discovery that the reality of the free energy inevitably clashes with "physical" con- 
straints arising from matrix- mo del/topological field theory Ward identities Q. The PI 
transcendent was subsequently demoted to the status of a generating functional for the 
string perturbation series. This fate has been shared by all solutions to the string equa- 
tions for gravity coupled to unitary c < 1 conformal matter. Many attempts to circumvent 
these difficulties have been made. Those which are physically well-motivated have proven 
to be mathematically inconsistent 0. Others, while mathematically consistent lack any 
cogent physical rationale. We are thus left with the unsatisfying situation of having no 
nonperturbative unitary c < 1 theory of gravity. 

There is a two- fold origin of the above dilemma. First, there is no simple spacetime 
interpretation of the c < 1 models. Hence questions of, e.g., whether the specific heat 
should or should not have an imaginary part are difficult to resolve. Second and more 
importantly, there is no physical principle which isolates a reasonable parametrization of 
physically acceptable possibilities. 

The first difficulty is not essential thanks to the c = 1 model of 2D gravity [Ej.0 



In the c = 1 model there is a simple spacetime interpretation |K|: strings move in two 
target space dimensions. The nongauge field theoretic degrees of freedom of the string are 
described by a single massless boson field - the "massless tachyon" which is related to the 



eigenvalue density field of collective field theory [|lTi |12l. The vertex operator calculations 
of the c = 1 matrix model are Euclidean continuations of the scattering amplitudes of the 
massless boson in a half-space, i.e., the "wall" S'-matrix amplitudes in the terminology of 



Polchinski 13 



With this clear physical interpretation there is an equally clear physical constraint 
on the theory: the ^-matrix must be unitary. That this is true perturbatively might be 



1 For reviews see 0-^). 
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expected to follow automatically either from the point of view of collective field theory 
or from string perturbation theory. However the question is somewhat more subtle for 
nonperturbative definitions of the theory. 

In the present paper we study nonperturbative unitarity of the c = 1 S-matrix. We 
will show that the situation at c = 1 is the reverse of that at c < 1: while nonperturbative 
unitarity may be used to rule out some theories there is a plethora of nonperturbatively 
unitary c = 1 theories. Thus, the second and more important difficulty alluded to above 
will be exacerbated. Indeed, our construction applies to a wide class of matrix model 
potentials V(X). These may be divided into two classes by the asymptotic properties of 
V(X) as the eigenvalue variable A — > — oo. In theories of type I V(X) — > +oo rapidly for 
A — > — oo and V(X) — > —A 2 for A — > +oo. Thus the theory is effectively defined on a 
semi-infinite line. The canonical example of such theories is defined by V(X) oc —A 2 on the 
interval [A, oo) together with an infinite wall at A = A. We will sometimes specialize to 
the case A = where formulas simplify. Although the potential is not analytic at A = A 
we expect that for a smoothed out version of the wall the main results will be unchanged. 
Theories of type II are defined by a smooth potential V(X) — > —A 2 on both ends of the 
real A axis. Thus, in perturbation theory there appear to be two disconnected worlds. 

In outline, the paper is organized as follows. In section two we give a definition of the 
S'-matrix in terms of large spacetime asymptotics of correlators of the eigenvalue density 
operator of the matrix model. We begin with an integral representation for the correlators 
derived in [IJJ] and continue the resulting asymptotics to Minkowski space to obtain the 
S'-matrix. The procedure of this section is an important technical advance over previous 



calculations. For example, in [14], the two-, three-, and four-tachyon scattering amplitudes 
were obtained to all orders of perturbation theory using the small length asymptotics of 
macroscopic loops. However, at the time |14j] was written, extension of the results to more 



general amplitudes appeared hopeless. Moreover, while the results of [|14| are valid to all 
orders of perturbation theory the nonperturbative foundations of these formulae are shaky. 
Another advantage of the procedure of section two is that the nonperturbative formulation 
of the theory is unambiguous. 

In section three we formulate the result of applying the procedure of section two in 
terms of some simple graphical rules. These rules lead to a relatively simple and compact 
form for the n m amplitudes (eq. (3.6) below). Some particularly simple cases, for 
example the nonperturbative 1 — > n amplitudes are written out explicitly (eq. (3.11)). 



The graphical rules also provide considerable insight into the fundamental nature of c = 1 
scattering, making clear how particle production is possible in a theory of free fermions. 

In section four we apply the lessons learned from the graphical formalism of section 
three to the problem of nonperturbative unitarity. Both theories of type I and theory II 
are unitary to all orders of perturbation theory. In section (4.1) we prove that theories 
of type I are nonperturbatively unitary. In section (4.2) we remark that theory II is not 
unitary, essentially because the theory does not take into account fermionic soliton sectors. 
Our remark will strike many readers as trivial. Nevertheless, we think it is important and 
deserves emphasis. The unitarity proof for theory I suggests various interesting general- 
izations of c = 1 scattering and leads to a characterization of a large class of acceptable 
nonperturbative definitions of the c = 1 .S-matrix. 

The exact formulae for c = 1 scattering amplitudes allow us to investigate in some 
detail the analytic structure of the S'-matrix in section five. We find several interesting 
singularities of the analytically-continued ^-matrix elements and interpret these in terms 
of metastable bound states of matrix model fermions. We emphasize the dependence on 
the parameter A. This parameter does not appear in perturbative amplitudes but has an 
important influence on the nonperturbative answers. 

In section six we make some remarks on the worldsheet/Liouville interpretation of 
our results. Most importantly, we show that the vertex-operator-motivated prescription 
used in |T^] to extract correlators from small loop-length asymptotics of macroscopic loop 
amplitudes is equivalent (to all orders of perturbation theory) to the procedure of section 
two. 

In section seven we make a few remarks on the background-dependence of the S- 
matrix, focusing on the background dependence of the singularity structure uncovered 
in section five. The dependence on tachyonic perturbations of the background will be 
addressed in a separate publication. In section eight we study some aspects of the high 
energy behavior of our scattering amplitudes. The amplitudes exhibit many nontrivial 
features. While this is hardly surprising from the close relation of the present S'-matrix to 
that of potential scattering, the interpretation of these features might prove instructive for 
string theory. In the conclusion we present some remarks on directions for future research, 
and we summarize some technical details in several appendices. 



2. Definition of the S'-matrix 



2.1. Eigenvalue Density Correlation Functions 

The Euclidean Green's functions of the eigenvalue density p = ftip(\,x), where x is 
the "time" dimension of the c = 1 matrix model are defined by 



GWlidean^l, A 1} . . . , X n , A n ) = (fJ,\T ^ 1p (x i , X X ) ■ • • ft 1p(x n , A„) J |/i)< 

GWiidean(<Zi, Ai, . . . q n , A n ) = J Y[ dxie iqiX% G{x x , Ai, . . . x n , A n ) . 



(2.1) 



Since the fermions are noninteracting these Green's functions may be written in terms of 
the the Euclidean fermion propagator: 



S E (x 1 ,\ 1 ;x 2 ,\ 2 ) = e-» A * ^e-^IfoAi.Aa) 

J-oo ^ 



(2.2) 



where / is the resolvent for the upside down oscillator Hamiltonian H = hp 2 — h\ 2 . In 



particular, for q > 

I(g,Ai,Aa) = (J(-g,Ai,A 2 ))* = (Ai 



1 



|A 2 ) 



H — p — iq 

Thus we obtain the integral representation for the eigenvalue correlators |I3 



(2.3) 



CEuclidean(q , i, Ai) = - / II T^" 6 "^^ /J JJ S E ( x a(i) j A CT (i)5 , A CT (j +1 ) 



71 J ~°° o k=l 

where Ql = q + q a{1) -\ q a{k) . 

2.2. Relation to Collective Field Theory 



(2.4) 



i). 



We now relate the exact eigenvalue correlators (2T1) to the correlation functions of the 
Das-Jevicki collective field formulation, as explained in pT^,|T3|,p~5|,|T6|| .§ 



Of course, the connection had also been extensively discussed in Jl2,16,13] before [15]. Nev- 



ertheless, the method of appendix A in [15] leads to the simple calculation of this paper, which 
we believe is new. 
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First, let us find the large A behavior of /. We work on half of the A axis, change 
variables to A = 2 v //Icoshr and identify dp = p — (p) = d\x as usual. If we hold r 
fixed then asymptotically at large p the integral I has the form of a direct and reflected 
propagator [ |T5|| : 



I{q,Ti,T 2 ) ~ 



\J Ap sinh n sinh r 2 



e -g|r 1 -r 2 | e i M |G(r 1 )-G(r 2 )| jD ^ ) 



+ie -9(ri+r 2 ) eJM (G(r 1 )+G(r 2 )) jR ^ j ^ ^ 



(2.5) 

Here G(r) = r — \ sinh 2r is the WKB phase factor of a fermion wavefunction, and R = 
1 + + J 2 -C'( 2 ) + • • •. The C« are real, polynomials in the g^, and rational functions 

of Wi = e Ti . Finally, D(q, t\, r 2 ) = R(q, T\, -r 2 ) for n > r 2 . 

The correlation functions of the Das-Jevicki collective field theory are obtained by 
inserting ( |2.5| ) into ( |2.4[ ), expanding the product, and keeping only those terms for which 
the nonperturbatively oscillating factors e ±lMG ( T ) cancel. This defines a new set of Green's 
functions G(qi,Ti;p) (at least as asymptotic expansions in l/p) which are the Euclidean 
continuation of the Green's functions of the Das-Jevicki theory to all orders of perturbation 
theory: 

(2.6) 



G r Euclidean(T'i 5 qi', A 4 ) — J][ 



1 



-<oin^rf>< 



2 yfp, sinh Tj 

Although in principle this procedure gives a straightforward method for calculating 
the off-shell Green's functions of d T x, in practice it becomes cumbersome already at one- 
loop [0. Nevertheless, as we will see, it is an extremely useful observation for extracting 
the ,S-niatrix. 



2.3. Large A Asymptotics 

Our strategy for extracting S'-matrix elements (in Minkowski space) and (microscopic) 
vertex operator correlators (in Euclidean space) will be to extract the coefficients of appro- 
priate terms in the large A^ asymptotic expression of G. From the collective-field-theory 
we may invoke a coordinate space version of the LSZ reduction procedure to identify the 
S'-matrix as the coefficient of the term multiplying appropriate incoming and outgoing 



wavefunctions of r |Tr8|,|19 

In the r 
Since 



00 limit the complicated rational expressions in (|2.5|) simplify drastically. 
) is just the Green's function for H it may itself be expressed in terms of parabolic 



5 



cylinder functions. From this representation one obtains the asymptotics for Aj — > ±00. 
The calculation is outlined in appendix A. For Aj — > +00 we obtain 



^(?,Al,A 2 ; 



— i 



vAiAi 



e -g|ri-r2| e iM|G(ri)-G(r 2 ) 



e ^(G(r 1 )+G(T 2 )) e -g(r 1 +T 2 ) 



(l + O(e" T 0) ?>0 



/(?,Al, A 2 ; 



vAiAi 



(2.7) 



e q\ri-T2\ e -i/J.\G(Ti)-G(T 2 ) 



_|_^ )* e -V(G(Ti)+G(r 3 )) e g(Ti+Ta) 



(1 + O(e- T 0) ?<0 



The "Euclidean fermion reflection matrix" ir^ will play a key role in what follows. In 
theories of type I, defined on A G [A, 00) we find 



9 P V 1 - ie- w ('*+*l«l) V r(i |g|) 

for A = 0. For nonzero A the result is given in equation (A. 7). 
In theory II we have 



(2.1 



R q — 



(2.9) 



The factors R in ( |2.8| ) and (|2.9| ) agree to all orders of perturbation theory 

1 



R 



11 



-Ri 



(2.10) 



9 1 _|_ ie -K{v+i\q\) 1 

and have a perturbative expansion of the form R q ~ 1 + 0(1/ /a). In theory II there are 
two worlds and we must also calculate the "transmission" probabilities obtained from the 
asymptotics of I for Ai = 2 v /^coshri — > +00 and A 2 = — 2 v /^cos1it 2 — > —00. In this case 
the "direct propagator" in (|2.7|) does not appear and the transmitted amplitude is 



(2.11) 



We may now derive the large A asymptotics of the collective field theory Green's 



function G. Substitute (|2.7|) into Q2.4Q and expand, keeping only terms for which factors 
of e lfJ,G ( T ^ cancel. For large r the two-point function behaves like: 



^Euclidean = <5(<?i +<? 2 ) 



1 



4/i sinh n sinh r 2 



\ q \ e -\i\hi-T 2 \ + K 2 (qi, - qi ; n)e-^ iTl+T2) 



(2.12) 
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giving a direct and reflected propagator for the string-theoretic massless tachyon. Specifi- 
cally, for q > 

TZ 2 (q, -q; //) = F dxR x R* q _ x (2. 13) 

Jo 

For three or more operators the resulting expression for G will have large r» asymptotics: 



e -\Qi\ri 



1 + C(e- Tl ) 



(2.14) 



Equation ( |2.14j ) may be taken as the definition of the functions lZ n . In theory II we define 



in an analogous way the functions lZ n (ei, q{) where = ± for \ — > ±oo. 

2.4- Continuation to Minkowski Space 

The analytic continuation to Minkowski space is most clearly understood by starting 
with the Minkowski space formulation of the Das-Jevicki theory, continuing that theory to 
Euclidean space and comparing with the Euclidean free fermion amplitudes. By studying 
the behavior of propagators we can obtain the following analytic continuation prescription. 
Consider the function H n (qi, . . . , q^\ q[ . . . q[) in some kinematic region where qi > and 
q[ < 0. We replace qi — > — iu>i and q[ — ► iu^, or, briefly, |g| — ► — iiv. Here and hereafter ui 
represents a real positive number, the energy of a massless tachyon. 

Explicitly, in theory I with A = the reflection factor becomes a pure (energy- 
dependent) phase: 



q p V 1 -ie-*o*«) y r(| + z(^ + ^)) p 

i? g ^i^lo^^e-i®^-") q<0. 

(2.15) 

More generally the analytic continuation will define a phase e * e (M+^,^) The factors of 
i e ifi,logfi an( ^ ^ioj ^ no ^. a ff ec ^ final amplitudes. The former cancels out of eq. (3.7) below 
while the latter may be absorbed as a phase redefinition of the states. 

In theory II R does not become a pure phase but has absolute magnitude (1 + 
e -27r(^+a))^-i/2_ Similarly T has magnitude (1 + ^(11+^)^-1/2 _ p Qr ^ } ar g e anc j p OS itive 

any S'-matrix element involving k bosons traversing the barrier is exponentially suppressed 
and of order e - 2nk v_ 
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2. 5. Definition of the S -Matrix 

In theory I there is an incoming and outgoing Fock space of bosonic states defined by 
oscillators a(u) and normalized by 

[a(u), a\uj')] = uj5(uj - u') . (2.16) 

In the space of incoming states 7i in , \uj\ 1 ...Uk) represents a set of left-moving bosons 
approaching the wall from the left with wavefunction 

tik^iu,*) = -±=Y,Xl e ~^ ){u+Tz) ■ ( 2 - 17 ) 

* ■ a 

In the outgoing space {u\ represents a set of outgoing right-moving bosons with wavefunc- 
tion 

= ^Eri e "^ (i)(t: " r:) • ( 2 - 18 ) 

The wavefunctions ip R , tp L are continuations of the Euclidean wavefunctions Y\ i e~\ qi \ Ti+tqiXi . 
Thus, the analytic continuation of the previous subsection shows that negative and positive 
q correspond to incoming and outgoing particles respectively. 

Putting together ( |2.14| ) with ( |2.17| ) ( |2.18| ) we may write the large r asymptotics of the 
Minkowski space Green's function as: 

k I k I 

SMtafamBkifoT)"^ (2.19) 

i=l i=l i=l i=l 

■k + l 

where S c is obtained from -^=]Z n by analytic continuation. 
Equation (|2.19[ ) defines the connected ^-matrix element: 

(uilSlu't) = d(^^uj i — ' s } j uj' j )S c {oJi\oj' i ) + disconnected terms . (2.20) 



3. Calculation of the S'-matrix 

3.1. Diagrammatic Formalism 

The procedure of extracting terms for which e ±ZAtG ( T ) cancel from the product of the 
functions I can be given a diagrammatic interpretation which facilitates calculation and 
which will be crucial to our proof that theory I is unitary.! 

3 The formalism is closely related to old-fashioned non-covariant perturbation theory for the 
fermions. Several special features of our problem lead to ad hoc rules which necessitate our detailed 
discussion below. 
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As in a Feynman diagram there is a vertex in the (x, r) half-space corresponding to 
each operator i(j*ip(x,T). While the final result will of course be independent of the order 
in which the Tj are increased to infinity, in intermediate steps we will choose some order 
and locate the vertices accordingly. Points are connected by line segments, representing 
the integral /, to form a one-loop graph. Since the expression for / in (2"7F) has two terms 



and we have both direct and reflected propagators as in fig. |l|. Each line segment carries a 
momentum and an arrow. Note that in fig. [I] the reflected propagator, which we call simply 
a "bounce," is composed of two segments with opposite arrows and momenta. These line 
segments are joined to form a one-loop graph according to the following rules: 

RH1. Lines with positive (negative) momenta slope upwards to the right (left). 
RH2. At any vertex arrows are conserved and momentum is conserved as time 
flows upwards. In particular momentum qi is inserted at the vertex as in 

fig. a 

RH3. Outgoing vertices at (x ou t,r ou t) all have later times than incoming vertices 

(■^ini ^"ire)' ^out ^ 



J in • 



Since diagrams drawn according to these rules correspond to real processes taking 
place in spacetime they will be called real histories. Some examples are shown in fig. |3| 
and fig. |j. There is a finite number (less than n!) of ways of connecting the dots. Since 
^2 qi = there is an overall undetermined loop momentum q, however the constraints 
RH1-RH3 restrict q to lie in a finite interval. 

To each real history we associate an amplitude, easily derived from (|2.7|) and illustrated 
in fig. [5| and fig. |6|. The total amplitude is simply obtained by summing over real histories 
to produce a formula which reads, schematically, 

RH bounces 

In theory II we take two copies of (x, r) space. Transmission amplitudes are computed 
by connecting vertices in one half space to another. Each transmission line is weighted with 
a factor of T q . Note in particular that since two fermions make a boson every transmission 
of a boson has a nonperturbative suppression of e _27I " M . 

9 



3.2. Explicit Formula 

We now write the equation ( |3.1| ) more precisely. Each real history has the structure 
of a series of "brackets" fig. [7] separated by bounces since (by considering momentum 
conservation) direct propagators can only connect incoming to incoming and outgoing to 
outgoing vertices. Each bracket must contain at least one insertion of momentum, so in 
ann-»m amplitude there can be up to Min{n, m} incoming and outgoing brackets. The 
number of incoming brackets equals the number of outgoing brackets. As we traverse the 
loop, we encounter an increasing set of momentum insertions. The factors for each bounce 
only depend on the net momenta flowing through the bracket, but there is a combinatoric 
factor counting the number of ways of forming a bracket out of a given set T of momenta. 
Consider for example the case where T is a set of positive momenta. The only constraints 
on the formation of a bracket as in fig. |7] is that, if q* corresponds to the vertex r* with the 
largest r in T then the direct propagators immediately before and after the insertion of q* 
must have positive and negative momenta, respectively. Let Tj be the momenta inserted 
before the insertion of q* and T 2 the momenta inserted after. The sum over real histories 
involves a sum with weights: 

E 0(q(T 2 ) + q*- Q)0(Q - q(T 2 ))(-l)W 
T 1 UT 2 U{q*}=T 

E (W - Q(T 2 )) - 0(Q - q(T 2 ) - <f )) (-1)^ 

TiUT 2 U{g*}=T 

= E (-^(Q - q(t 2 )) 

TinT 2 =T 

= - E (-i) |T2l %(^)-Q) 

TiUT 2 =T 

= f+{T t Q) 

(3.2) 

where the notation q(S) means the sum of momenta in a set S and in going from the 
second to the third expression we have noticed that the two terms in the sum correspond 
to subsets which do and do not include q* . Note in particular that the final weighting 
factor makes no reference to q* and is therefore independent of r-ordering. Similarly, for 
a set T of negative momenta as in fig. [7] we have a weighting factor 

f-(T,Q)= E %(Ti)-Q)(-l) |T21 

TlUT - T (3.3) 
= - E W-9CTi))(-l) |Ta| • 

TiUT 2 =T 
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Note that the form of / differs for positive and negative momenta. For convenience we set 

/±(0,Q) = o. 

We are now in a position to write equation (|3.1| ) precisely. Let S = {q±, . . . , q n } be the 
set of momenta, the set of positive and negative momenta, and define an admissible 
filtration (AF) of order k to be a tower of subsets of momenta: 

= F C F x C F 2 c • • • C F 2k = S (3.4) 

such that 

Fj - Fa i CS+ j = (mod 2) 

(3.5) 

Fj - Fj_i C S~ j = 1 (mod 2) . 
Then we have the (Euclidean space) "filtration formula" 

777< j n k 

^ = r E/:E / ^II/(-iy( F i- F i-i.« + 9( F i))n^(^-o^(i^) ( 3 - 6 ) 

k=l AF k j=l 3=1 

where m = Mm{\S + \, \S~\}. The combinatoric factors f± impose kinematic constraints 
which render the Q-integral finite. 

Examining the reflection factors of (|2.8|) , (|2.9|) we see that the fx dependence of R q is 
only through the combination fi + i\q\, except for the (noncancelling) prefactor jU - ' 9 '. This 
observation allows us to compute the derivative (/i2 S ^7V) by converting it to a Q 
derivative. This has the effect of killing the loop integration, extracting a boundary term 
wherever the argument of one of the theta functions vanishes. Thus we can write a formula 

■fL^ M Hn) =i n+1 ^S |9il E ± II R Q R *Q ( 3 - 7 ) 

" ^ ' RH' bounces 

where the real histories allowed in (13.71) have the added condition 



RH4. Exactly one direct or reflected propagator carries zero-momentum. Direct 
propagators with zero momentum have a factor +1. 
The following useful properties are immediate consequences of these expressions. The 
functions TZ n (qi, . . . , q n ; fi) are real and totally symmetric functions of the qi and are 
defined on the plane = in M n . If S is any subset of momenta then the functions 
are continuous but have discontinuous derivatives across the hyperplanes H$ defined by 
q(S) = J2 q eS 9 = 0. On the complements of these hyperplanes the lZ n are analytic. 
Indeed, we can Taylor-expand for small qi and the Taylor-coefficients are polynomials 
in polygamma functions (see e.g. appendix C), although the coefficients in the Taylor 
expansion change across Fig. Similarly, in the asymptotic expansion of lZ n for \i — ► +oo 
the coefficients of 1//U n are polynomials in the qi, although the polynomial changes across 
the hypersurfaces Hs- The analytic structure of lZ n will be discussed further in section 5 
below. 
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3.3. A Low-Energy Theorem 

Suppose any momentum is taken to zero. Physically the vertex operator inserting 
becomes proportional to the cosmological constant and the resulting limit should be 
expressed as a derivative with respect to \i of an n — 1 point function. This can be proven 
by noting that if G Fj — Fj_\ then 



f±(F j -F j _ li Q) -> -p* ^f±{F' 3 -Fj-uQ) (3.8) 



where F'- = Fj — {p*}. Hence combining this remark with ( |3.7|) 



^l^mnn^p.— i^^^TZn-! ) (3.9) 



as any momentum p* goes to zero. 



3.4- Special Cases 



The general formula ( |3.6| ) is somewhat awkward to work with. We give the explicit 
formula in some special cases. 

1. If S + = {qi, . . . q n -l}, Qn < 0. Then we have: 

§-(^K)=i n+l £ (- 1 )' 52l ^ 1 )^(5 2 ) • (3-10) 

Correspondingly, in theory I the S'-matrix element for an incoming state |u>) to produce 
(ui, . . .,uJ n -i\ is given by 

•71+1 ru(S 2 ) 

V^Y- Sl -- Jo 



us 2 =s+ 



This allows us to compute the probabilities for particle production when an incoming 
tachyon impinges upon the wall. A discussion of this effect follows in section 8. 

We also note that the content of ( |3.10| ) and ( |3.11|) can be expressed in terms of linear 



Ward identities on the amplitudes, which could be interesting in Liouville theory. 

2. We can use our rules to write a relatively simple result for the S- matrix for scat- 
tering ui + U2 — > E w i- This consists of two terms with either 2 or 4 bounces. Let 
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S + = . . .q n -2}, and S = {<? n _i,<? n }, then: 



E (-1) |T21 / " (l - 0(Q + 9 n-i) - 0(Q + q n ))R* Q R Q+q(s -)dQ 



TiUT 2 =S+ 



+1 E 



SiUS 2 =S+ 



q(S 2 ) 



f+(S 2 , Q)R q ^s 2 )-Q R *Q 



f+(Si, Q + q n + q(Si))R qn+ QR* Q+qn+q{Si) 



+ f+{Si,Q + q n -i + q( s i))Rq n -i+QR*Q+ qn _ 1+q (s 1 



(3.12) 



Specializing further, consider 2^2 scattering. Let n = 4 above, and take qi+qs > 0, 
qi + q± > 0. We find: 



-4//2EI 



1 Jm 



Rqi+q2R() Rq±R 



f/2 



D D* p p* p p* p p* 

- n '92- n -(j2+'?3 - n '<?4- n -0 - rL g2- rL g 2 +g4 ^-"-O 



(3.13) 



The previously published result for 2 — > 2 scattering in |TJ[] involved an infinite sum of 
gamma functions and was rather unwieldy. This illustrates nicely that the formulae of 
this paper are a substantial improvement upon those given in ||14|1 . (And also that there 
are some remarkable identities on gamma functions, this particular case is written out in 
appendix B.) 



4. Unitarity of the S'-matrix 

4-1. Theories of Type I are Unitary 

The essential idea of the unitarity proof is to regard the time evolution of the real 
histories of the previous section as a composition of three maps: fermionization, free 
fermion scattering, and bosonization: o Sff o as in fig. |[ 

The most complicated map is the bosonization map, although this is well-known. We 
will let a\v),a{y) be delta function normalized fermion creation operators acting on the 
Fermi sea: 

a(e,+)|/x> = at(// + e)|A*> = £>0 
a(£,-)|/i) = a ( M - e)|/x> = £ >0 . 



(4.1) 
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An operator in the fermionic theory is normal ordered if it is a sum of terms of the form 
Hi a (Ci? e i) with e^i < 0. Under bosonization a one-particle incoming state is mapped 
according to 

/oo 

(4-2) 

= / d£a{fj, + Z)(J{fj, + t-u>)\n) . 
Jo 

Note that if incoming states are normalized according to (co\u') = u5(u — uj') then the 
map ►/ is an isometry. This isometry may be extended to the entire Fock space by 

/oo n 
Y[d& Y[a(n + £i)a^ (fi + & -Ui)\fi) . (4.3) 
-°° i i 

Note that the operator acting on \fj) on the RHS is not normal-ordered and hence the 
state does not have definite particle-hole number. We may now imagine normal-ordering 
the operators in Q4.3|) . To each normal-ordered fermionic monomial we can associate a 
diagram of the type drawn in our real histories. In particular, the state Y\ i a(£j, €j)|/i) 
created by a normal-ordered operator is associated with a set of incoming bounce-lines 
carrying negative momentum and having upward (downward) - pointing arrows for 
e = ±, respectively. By an inductive argument one can show that the normal-ordering 
process corresponds exactly to forming all possible incoming bracket structures discussed 
in subsection 3.2 above. The details of this argument appear in appendix D. 

Thus the first third of the set of real histories corresponds in a precise sense to the 
fermionization of the incoming bosons. The second stage corresponds to free fermion 
scattering. In theory I this scattering is diagonal in the momentum basis Yl i a(£i, ej)|/x) 
and merely reverses the signs of all ej. In particular, 

n m n 

(n\l[a(Z h e i )S ff l[a(Z' j ,e' j )\n) = 8 n , m £ ]J Sfa + Cu))^.^®^ . (4.4) 

i=l j=l o-6£„j = l 

The final third of the map is just the inverse of the isometry £&—►/, since our graphical 
rules implement this inverse (simply consider reading the diagram from top to bottom). 
We have written the bosonic S'-matrix as a composition of two isometries and a unitary 
free-fermion ^-matrix; it is therefore nonperturbatively unitary. 

Remarks: 
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1. The above discussion allows us to characterize a large set of nonperturbatively 
unitary c = 1 S'-matrices. Note that at no point were any special properties of the function 
@(x) employed: the unitarity equations may be checked without resort to any identities 
on the Bernoulli numbers. The specific form of O only enters in the comparison with 
string perturbation theory. Thus we may construct a nonperturbatively unitary string S- 
matrix (in two dimensions!) using any real-valued function Q(x) which has an asymptotic 
expansion 



Q(x) ~ argT(^ 



IX 



2 

oo 



~ -xlogx + X + J2 hf^iV 1 - 2- (2 ^ 1) )x- (2n - 1) (4 ' 5) 
' ZniZn — 1 

n=l v ' 

for x — > +oo. For example, putting the infinite wall at a position A = A leads to the 
reflection coefficient given in eq. (A. 7) below. For any A this coefficient defines a distinct 
S'-matrix, corresponding to a Q(x,A) satisfying ( |4.5|) . Indeed, by "inverse scattering" 
we expect one can assign a matrix model potential V(X; 0) to any function O satisfying 
( |j~5|) . It would be interesting to describe the set of potentials which are characterized by 
the requirement ( [4.5| ), since, according to the philosophy of section 7 below, this ambigu- 
ity corresponds to the existence of differences of background geometries detectable only 
nonperturbatively. 

2. The fermionic formulation partially clarifies the question of the symmetry 
which the system is known to respect [p0| - ^0i[T5[| . This is not manifest as a symmetry of the 
S'-matrix we have computed. From the point of view of the free fermion system, the Cartan 
subalgebra of generated by 2s, ° is a manifest symmetry of ( |4.4[ ), corresponding to the 
infinite set of conserved charges Q2 S = z\2 P 2s constructed from the fermion momenta (note 
that due to the reflection off the wall, only even powers are conserved). The realization 
of the rest of the symmetry algebra is not obvious. We now see why the symmetry is 
so deeply hidden in the bosonic formalism - the conserved charges do not commute with 
tachyon number. 

3. The above formalism also makes clear how it is that one can have particle produc- 
tion in a "free" theory. The essential point is that the fermions scatter with an energy- 
dependent phase e l ®( E ' v \ Thus a particle-hole pair corresponding to a single boson scat- 
ters into a "dispersed" particle-hole pair with overlap on states of arbitrary boson number. 
This may be regarded as a nonperturbative generalization of the dispersion of tachyon 
wavepackets described in [|13| . 



15 



4. The above discussion suggests an obvious yet interesting generalization of the 
present S'-matrix. We may imagine adding "flavor" quantum numbers to our bosons and 
fermions and may further replace the crossings of fermion lines "in the bulk" by any 
factorizable S'-matrix compatible with the bounce factors. The compatibility conditions 
have been investigated in []31],[32] . The above remarks also provide a general mechanism by 
which one could use factorizable S'-matrices to write nontrivial S'-matrices with particle 
production. 

4-2. Theory II is not Unitary 

The nonunitarity of theory II is a simple consequence of the fact that an incoming 
particle-hole pair can dissociate as in fig. |] into a leftmoving and rightmoving fermion. 
One fermion can be reflected back into its original half space while the other can tunnel 
through into the adjoining half space. The Hilbert spaces for left- and right- movers are 
then the one-soliton sectors, which have no analogue in the Das-Jevicki theory.0 

More formally, we may consider the unitarity equations for 1 — ► 2 scattering on the 
same side of the potential. The contributions of transmitted amplitudes to the equation 
must be of order e _47I " M . However, the reflected amplitudes fail to satisfy the unitarity 
equations at order e _27rAi , as we now show. Expanding in small energies u>i we have the 
S'-matrix elements 

S(u, +) = ^ — u H 



S(u, +\uJl, +; L02, +) = —77 7"2 ^1^2 



1 7T 



S(ui,+;u 2 , +) = — — UJ1IV2 



TV 1 

2 cosh 2 t\[x 

+ iu 1 w 2 (u 1 +u 2 )( - - \ 2l ,A \ + 77 73 i>o) +■■■ (4.6) 

1 + e ^ 2 cosh ix n 

7T 1 

2 cosh 2 -Ki± 

+ iuiu 2 {ui +u 2 )( - -32^ ^1 + 77 — A V'o) +••• 

1 + e ^ 2 cosh irfi 

where ipo = Re(V'(^ — i/J>)) and ipi = The first term is order (9(e _27rAt ) while the 

second is order 1. Checking the unitarity equation (u, +\SS*\u>x, +; u> 2 , +) = in a small 



4 



The Das-Jevicki theory does have soliton sectors corresponding to one-eigenvalue instantons 
[33]. It is possible that taking these into account will restore unitarity within the context of the 



collective field theory. 
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energy expansion we see that the 0(uj 3 ) terms, which are perturbative, cancel but the 
nonperturbative 0(uj 2 ) terms cannot cancel. Moreover the contributions of the transmitted 
amplitudes to the unitarity equations are 0(e~ 4nfl ). Hence unitarity is violated at order 
e -27T/x ag eX p ec t ec [ f r om the simple physical picture. 

The violation of unitarity in the massless tachyon theory is extremely suggestive in 
view of the conjectured connection of the c = 1 matrix model to black hole physics [j34 -j37|. 
Unfortunately, that relation must be further clarified before we can confidently apply the 
above results to the black hole problem. 



5. Analytic Structure of the S-matrix 

5.1. General Remarks 

In standard field theory the analyticity properties of an S-matrix are related to im- 
portant physical properties of the theoryi For example, causality implies analyticity in 
suitable domains, and the existence of unstable particles and thresholds imply the exis- 
tence of poles and cuts in analytically-continued S-matrix elements. Therefore, in this 
section we examine the analytic properties of the c = 1 S-matrix hoping to understand 
better the nature of the physics of 2D strings. 

Many aspects of our problem are different from the more familiar examples of rel- 
ativistic field theories in Minkowskian spacetimes. A first difference is that in stan- 
dard S'-matrix theory [|38| S'-matrix elements are considered as analytic functions of 
Sijk... = ipi ±Pj + • • - ) 2 - Here we will analytically continue in the energies u>i, which 
are, roughly speaking, y/Sijk...- A second (related) difference is that in our problem we 
only have time translation and time reversal invariance so many restrictions of Lorentz in- 
variance are lost. A third difference is that the particles in our case are massless, whereas 
in more realistic examples massless particles like photons and gravitons, which interact 
with one another in asymptotic regions of spacetime, are usually excluded from discussion. 

The third point above leads to a rather different analytic structure of S-matrix am- 
plitudes than that normally encountered in massive relativistic field theory. Consider, for 
example, the 2 — ^ 2 scattering amplitude in a scalar field theory with particles of mass m. 



In the complex s = (pi + P2) plane the amplitude has the structure indicated in fig. [10 



5 We would like to thank R. Shankar and A. Zamolodchikov for very useful discussions relevant 
to this section. 
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There are elastic threshold branch points at 0,4m 2 . On the physical sheet there are no 
other singularities, but if we analytically continue to the second sheet we may discover 
resonance poles at some complex mass-squared s = fj 2 as indicated in fig. [T(]. Moreover, 
the presence of these resonance poles implies the existence of further cuts, for example the 
one beginning at s = (m + n) 2 indicated in fig. [H]. Consider the limit of fig. PH| as m — > 0. 
In this limit the physical sheet separates into two disjoint half-spaces. There should be no 
singularities on the physical sheet and the physical S matrix is defined by two different 
analytic functions defined on these two regions. If we analytically continue, say, S + 
from the upper half plane into the lower half plane we will encounter resonance branch 
cuts located at the position formerly occupied by the resonance poles. 

In the c = 1 iS-matrix, the cut structure analogous to the separation of S + from 
S~ above has already been investigated at tree level in We will see below that 



nonperturbatively some new singularities arise corresponding to the presence of resonances. 

5.2. Analytic Properties of Bounce Factors 

The analytic properties of the S'-matrix follow from those of the bounce factors. For 
theories of type I with A = the bounce factor is 



«-■'-' \i r(i+ix) (5 _ 1) 



:y|e-/ 4 cos(f(I + zx))r(I 



IX) 



From the first equality it is obviously a phase for x real. From the second we see that it can 
be continued into the complex x-plane, where it has simple poles at x = zt = —i(2£ + §) 
for £ = 0, 1, 2 . . ., with residues pg = yfe t7r//4 (~l)^ (21+1)! • When A 7^ we have instead 

je(x,A) _ l + (e nx -i)X(x,A) iQ(x) 

l + ( e ™ + i)X(x,A) { } 
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In [39 1 D. Kutasov and Ph. Di Francesco observed an interesting analytic property of 



the tree-level 5-matrix. Namely, by restricting attention to "one-particle irreducible" elements 
one could work with analytic expressions (polynomials in the uJi) valid in all kinematic regimes. 
Unfortunately this property does not persist at higher loops (as one might guess from physical 
grounds). Using the explicit formulae above it is easy to check that the "one-particle irreducible" 
amplitudes defined in [p9| are not analytic in the momenta at one loop order. 
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where 

X(x, A) = l_r£ +i-)r£ - z-)^ 1 1 f 27 2 J . (5.3) 

The bounce factor in ( |5.2| ) defines a meromorphic function of x which has a sequence of 
poles at a; = zi(A), smoothly evolving from the poles at A = 0. We denote the residue at 
the pole by pi(A). 

The poles z^{A) in the bounce factors have a simple interpretation in terms of the 
free fermions. The denominator of (|5.2| ) vanishes for z such that the wavefunction (see 
appendix A for a definition) Xr( z i A) = 0- For Imz < the wavefunction x~r( z i ^) decays 
as A — ► +oo so the condition Xr{. z i A = A) = is exactly the condition to have a bound 
state energy z in the upside down oscillator potential with an infinite wall at A. The 
"energy" z is complex because this bound state is unstable. As A — > — oo there is a 
deeper and deeper well on the far side of the parabola and the resonant state with energy 
zg(A) becomes more and more long-lived, i.e., the imaginary parts of zg(A) go to zero. 
Indeed, the poles zi(A) approach the positive real axis and become dense there, eventually 
becoming a massless free fermion cut. (It is possible to obtain some analytic results on 
Zi(A) which confirm the above picture deduced on physical grounds.) 

5.3. Theory II and the Limit of Theories of Type I 

One may take the limit A — ► — oo to obtain the bounce factor of theory iT provided 
one works at Imz > 0: 

lim R I q (z,A) = R I q I (z) . (5.4) 
At Imz < this is not true and we have instead 

lim RfaA) = (l + e-^R^iz) . (5.5) 

A— > — oo 

The reason for the discrepancy is that as A — ► — oo a new cut for the fermion, the coales- 
cence of the poles zi(A), appears along the real axis. Indeed, for Imz = the A — > — oo 
limit of R 1 is ill-defined. 

The analytic structure of R n may be studied using the formula One finds a 

series of poles in the lower half plane similar to that found for R , leading to a similar 
analytic structure for the scattering amplitudes in the type-II theory. 

19 



5.4- General Properties of "Massless Tachyon" Amplitudes 

We may now continue the amplitudes in (|3l|) by choosing all U{ but one to be inde- 
pendent and continuing in those. Upon analytic continuation the integrals in ( |3.6j ) become 
contour integrals, and singularities arise when the endpoints of the integrals hit poles of 
the integrand or when poles of the integrand pinch the contour. Since we have written 
the integral as a sum of terms there will in general be several choices for how to continue 
the contours in the separate integrals. One can define a notion of a "physical sheet" by 
starting in a given kinematic region and analytically continuing into an appropriate half- 
space. Note that if a bounce line carries positive energy uj then the corresponding factor 
e ±!6((t±w) poles at uj = ^ffi — i(2£ + |). Therefore, from rule RH1, if we define the 
physical sheet to be included in the half-space Im(uj) > 0, there can be no singularities in 
the integral. On the other hand, it can happen that when we continue out of this half-space 
a pole in the integrand is forced to hit the endpoint of the contour integral. This results in 
a branch cut singularity. Moreover special configurations of uii can produce more compli- 
cated singularities, and continuation to further sheets can produce new singularities. Such 
singularities are indicative of new states and degrees of freedom in the system not directly 
seen in the physical S'-matrix. These points are best illustrated by explicit examples. 



5.5. Analytic Structure of S{u>\u>) 
In the 1 — > 1 amplitude 



S(uj\uj-A)= ! dxe ie{uj -^- x ' A)+ieifl+x ' A) (5.6) 
Jo 

the integrand has a set of fixed poles at x = —fx + ze(A) and set of moving (with uj) poles 
at x = uj — fx — Zi(A). The contour integral is therefore forced to hit these poles when 
uj = ±fx + Zi(A) and at these points the integral has a logarithmic singularity. We define 



the second sheet by choosing the branch cuts as in fig. 11. From the discussion of section 
5.2 we see that the branch cut singularities are naturally described by a simple physical 
picture: one fermion gets trapped in a metastable state. The other fermion bounces off the 
potential and can radiate and reabsorb massless tachyons, thus producing a cut singularity. 

It is interesting to proceed and continue to the third sheet by passing through one of 
the cuts of fig. [TI[ For simplicity we consider continuing through a cut at u = —fx + zi(A), 
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although a similar story holds for the cuts at ui = \x + zg(A). Comparing the contour 
integrals defined in the x plane in fig. [12] we see these differ by 



A) = A) + p e (A)e i@ ^- z ^' A ^ (5.7) 

where the superscript refers to sheet number. The extra factor produces a new set of 
poles on the third sheet at u = zg(A) + Z£i{A). These poles may be interpreted as a 
metastable state obtained when a particle-hole pair resonates. Since the fermions are 
free the "energies" Z£ simply add. Since both are trapped, neither fermion can shake off 
massless tachyons so the singularity produces a pole and not a cut. Note that when A = 
these poles occur at u> = —in for n = 3, 5, . . .. 

5. 6. Analytic Structure of the Four-point Function 

Proceeding along similar lines one can investigate in detail, e.g., the four-point function 
S(u)i, W4; A) for u>2 < W3,o;4 < u±. We take uj\ to be dependent and continue from 

a real subspace of 1R 3 . One finds a rather complicated singularity structure. The most 
interesting singularities appear to be double-poles arising from the four-bounce terms. 
These correspond to processes where an incoming tachyon (say, \oj%)) becomes trapped in 
a resonant state u>2 = zg(A) + zg>(A) which later decays to (u^. More precisely, letting 
UJ2 = Z£ + Z£> + e\ and = zi + Z£> + £2 with the e's small we get double poles 



Pi Pi' 



C i@(u> 3 -Zt) j 1 c iQ(w 3 -z p ,) 



£162 ei(ei-e 2 ; 



(5.8) 



The residue of this double pole, which is essentially a bounce factor may be interpreted as 
an amplitude for scattering off an "excited background" corresponding to adding a single 
resonant tachyon to the fermi sea. 

5. 7. Lessons from the Analytic Structure of S 

The main lesson we may draw from these remarks is that it is not possible to neglect 
the fermionic degrees of freedom beyond perturbation theory. Even in theories of type 
I, where we have a nonperturbatively unitary ^-matrix without including soliton sectors 
or matrix model fermions, the existence of the fermions can be detected by studying 
the analytic properties of the S'-matrix to find resonance poles and cuts. Presumably the 
nature of the residues and discontinuities at these singularities would indicate the fermionic 
nature of the particles. It must be emphasized that these analytic singularities are not 
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mathematical artifacts but have real effects. If A is large and negative the resonant states 
we have discussed have long lifetimes. If, for example, these lifetimes exceed those of 
experimentalists measuring the S'-matrix at A = +00 they will find themselves puzzling 
over an apparent loss of unitarity. 

A second lesson is that the rich singularity structure of the higher point amplitudes 
indicates the existence of a correspondingly rich spectroscopy of excited, but unstable, 
backgrounds. These time-dependent backgrounds, resulting from changes in the Fermi 
sea, are the nonperturbative analogues of the tachyonic backgrounds studied by Polchinski 
in fll3T 



6. Worldsheet Interpretation 

6.1. Connection to Liouville Theory 

So far our discussion has emphasized the free-fermion formulation of the matrix model. 
Since the double-scaled matrix model is a sum over continuum surfaces we expect that the 
continuum amplitudes can also be described by the conformal field theory of a massless 
scalar X coupled to a c = 25 Liouville theory <fi with worldsheet action 

A = J \dXOX + d<j)d<j) + V2~R (2) + (6.1) 



As explained in jlof the Liouville and r coordinates are not the same. Recall that the 
objects calculated in a sum over continuous geometries on 2-surfaces with boundary are 
"macroscopic loop amplitudes" defined by fixing the boundary values of the two-metric 
e -y/2^ go ifoai the bounding circles C have lengths £ = § c e^/^ [[40|--fl4 . 



In [|i4] , |44|1 it was proposed that one extract the c = 1 correlators of tachyons by 
extracting the terms proportional to nonanalytic powers £' p ' in the small £ expansion of 
the macroscopic loop amplitudes. From the continuum 2D path integral point of view 



W(£, p) corresponds to an expansion of a sum of local operators. As discussed in |15j that 
sum of operators can be written as 

W in (£,p) = -T p —^n-\p\/H lpl (2^jl£) -J2^, P jf z ^^ r/2 Ir(2^I£) (6.2) 

r=l 

where T p is proportional to the tachyon vertex operator and B rjP are redundant operators 
for p ^ 2. 
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We now relate the prescription of this paper to that of [n,|3|. From the matrix- 
model representation of macroscopic loop operators W(£, q) we see that the amplitudes 
are related by an integral transform: 



/"OO 

(JJW(^, 9i ))c= / \\d\ l e- i ^([{p{\ l ,q i )) L 

J A 



(6.3) 



in theories of type 1.0 The £ — > asymptotic behavior of the integral ( |6.3|) is dominated by 
the A — > oo behavior. Changing variables to r and using the asymptotic expression ( |2.14| ) 
we see that 



([[w{e u qj) 



?i2y / /iCOSh T, 



n 



(6.4) 



For small z and g^Zwe have the asymptotics: 

/oo 
rfre -£ 2v 7IcoshT e -MT _ (£^)kl r (_| ? |) 



(6.5) 



plus terms regular in £. Hence (11^) = T^ioXt • ■ ■ ? 5nj A 4 )- To make a connection to 
standard vertex operator normalizations we compare with the vertex operator calculations 
found, e.g., in |3^] which show that 



r q = IML [ e iix/V2 p V2(i-±\ q \)4> = r (M) 



r(-ki) h 



r(-M) 



V a 



(6.6) 



Thus we finally have a prediction for all integrated vertex operator correlators in the c = 1 
theory: 



V r (kil) 



f 1 \ 2ff-2 + n 



„ 3g — 3+n n 

Em / < n 6s n^.> 

g > \^/ JM g , n x i=1 



tj,=i 

(6.7) 

where M. g , n is the moduli space of curves of genus g with n punctures. 

(Let us note parenthetically that the "wavefunction normalization" factors r p^ 9 |j- > in 



7Jfy have been the subject of much discussion [[|5}|4^,|39| . These factors are singular when 



7 In theories of type II we must be more careful since the eigenvalue distribution grows on both 



sides of the axis. The procedure advocated in [14] is to take a Fourier transform with £ = iz, z real, 
and analytically continue positive and negative frequency components separately. This procedure 
is correct to all orders of perturbation theory, but its status as a nonperturbative definition is 
unclear. Hence we limit our discussion to theories of type I. 
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g G Z and it has ben proposed that the external line divergences should be reinterpreted 
as the contributions of "intermediate" on-shell states (in some space-time field theory). 
Such an interpretation requires inclusion of the "special state" vertex operators obtained 
from dressing nontrivial Virasoro primaries in Fock spaces with charge pe2 with highest 



2 

Virasoro weight ^V. An alternative interpretation of these factors was proposed in [I5| 



based on the smoothness of macroscopic loops for q e 2. The apparent singularity arises 
since the extraction of nonanalytic powers of I is not well-defined for 5 6 Z. Indeed the 
integral ( |6.5[ ) can be written more precisely as 

/ d Te -2VMCOShT e -| g |T = J J2y/Jie) 

Ja sm7r|g| 

» 1 e A(m-n-\q\) (6-8) 

- V(-inv^r v — r7 — tt n 

z — ' ^— ' m!(n — my. m — n — o 

n>0 m=0 v ' ^' 

Clearly the left hand side is smooth as q becomes integral, but the decomposition of the 
right hand side breaks down. This was interpreted in |15[] as the fact that £> Sj9 while in 



general a redundant boundary operator becomes a bulk operator at special values of q. It 
thus appears to be consistent to use the vertex operator normalization T q , bearing in mind 
that for g 6 Z the operator is a linear combination of special state and tachyon vertex 
operators.) 

6.2. Continuations to Minkowski Space 

The c = 1 x Liouville system is a conformal field theory with two Euclidean signature 
bosons. The Minkowskian spacetime physics deduced from this Euclidean theory depends 
strongly on how we analytically continue. 

The analytic continuation described in section 2 above identifies X as Euclidean time. 
Thus we continue X — > it and |A;| — > — iu. If we formally continue the expression in ( |6.7| ) 
we obtain correlators of macroscopic state operators ||!J 

e ±iut{z,z) e (V2+iu>)<t>(z,z) 

These operators have positive Liouville energy and, at least semiclassically, create surfaces 
with large holes and singular metrics. Indeed the difficulties associated with the so-called 
"c = 1 barrier" were ascribed by N. Seiberg to the destruction of a smooth world sheet 
arising from inclusion of such operators in the action [55].§ It would therefore be very 



See [50 1 for a different opinion. 
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interesting to see if the formal analytic continuation of the RHS of ( |6.7|) can be given a 
more rigorous justification. 

The analytic continuation of X to Minkowskian time also leads to a simple spacetime 
interpretation of the Seiberg bound [f49l , |5"0| which states that Liouville exponentials e a ^ 
must satisfy a < \Q. The fact that only one root of the KPZ equation is allowed is simply 
the fact that for scattering on the right half-line, incoming particles must be leftmovers 
and outgoing particles must be rightmovers. Note in this connection that if we work 
in theory II, with another asymptotic universe then, nonperturbatively, there can exist 
"wrong branch" states, i.e., incoming rightmovers from the point of view of the righthand 
universe.i 

Finally the issue of correct normalization of tachyon vertex operators discussed above 
becomes academic in this continuation |§. The S'-matrix obtained by continuing lZ n is 
simply related by the phase redefinition of statevectors 

\Ul . . . UJ n ) — > — ; \UJl . . . U n ) 

w . n ( 6 - 10 ) 
/ i TT T ( luJ i) i i 

- LJ - T(-iuJi) 

to the S'-matrix obtained from vertex operator calculations. Hence, no probability ampli- 
tudes are changed and the theories are physically indistinguishable. 

Another choice of continuation which is sometimes adopted is to regard the Liouville 
coordinate as Euclidean time. This prescription "0 — > i(p" is rather less well-defined. 
Indeed serious objections to it have been raised in IfjQf . Roughly speaking time becomes 



effectively semi-infinite, the Liouville wall corresponding to something like a big bang or a 
big crunch. In the former case the amplitudes TZ n (qi, . . . , q n ) correspond to "amplitudes" 
for producing right-movers and leftmovers (corresponding to positive and negative q^) 
from the finite past. Since time is semi-infinite there is no obvious analog to the unitarity 
constraint. 



9 A related remark was made by N. Seiberg. He noted that the wavefunctions of [44] associated 
to the "wrong branch" could be obtained from continuations appropriate to the far side of the 
parabola. 
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6.3. Mathematical Applications 

The equality of the asymptotic expansions in ( |6.7|) identifies the coefficients of the 
large \x asymptotics of 1Z n (pi, ■ ■ ■ , p n ', A 4 ) with integrals over M. g , n of some natural densities 
provided by the Liouville theory. From our construction of 1Z the answers are simply 
expressed in terms of Bernoulli numbers, thus giving formulae reminiscent of [ pT||52| . These 
observations add more evidence to the oft-quoted conjecture that there is a topological field 
theory version of the c = 1 model. To make the connection precise one should (1) generalize 
our formulae to the self-dual radius and (2) understand the appropriate nontrivial k — ► oo 
limit of the twisted N = 2 minimal model. Step (1) might be accomplished using the 
results of |53|,[5^1, but step (2) appears to be more difficult. 



7. Background Dependence of the S-matrix 

The most interesting questions in 2D string theory center on the nature of strings in 
different backgrounds. It should be possible to perturb the conformal field theory ( |6.1|) to 
obtain more general nonlinear sigma models with two target space coordinates {X 1 <fi). In 
general, nearby theories are expected to have an action of the schematic form 

A + 5A = J ldXdX + d(j)d(/) + V2R (2) (j) + lte^ + J dpe(p)V p + ^t r , a K, a (7.1) 

where V rtS denotes some basis of "special state operators." 

As long as the background perturbations preserve conformal invariance and the asymp- 
totically flat nature of the target space geometry, the string S'-matrix should make sense, 
and we should therefore discuss the functions S({ui} — > {<^}; backgrounds) on the infinite- 
dimensional space of background geometries. These can be parametrized, in an infinitesi- 
mal neighborhood of (|6.1|) by a subspace spanned by the directions e(p),t rjS . Ultimately, we 
would like to know what principle determines the allowed backgrounds and what equations 
govern the background dependence of S. 

The background dependence on e(p) is an interesting question which we will address 
elsewhere. For the present we describe how S changes in an infinitesimal neighborhood 
around (|6.1|) by perturbing the parameters t r>s . The matrix- model identification 

V r ,s <"> Bs,q=r = J AV^VKA, x)d\dx + ■■■ (7.2) 

26 



for r G {s, s — 2, . . . , — s} suggests, by exponentiation, that perturbations of the cr-model 
action by the parameters t r ^ s correspond to perturbations of the matrix model potential 

Moreover, this reasoning suggests that changes in V 
which affect the nonperturbative S-matrix but not its perturbation expansion could be 
interpreted as "small" perturbations in the background geometry which only affect non- 
perturbative physics. The dependence on A discussed in section five above is an example 
of such a perturbation. 

The change in the S matrix due to a change in background corresponding to 8V(X) 
is easily computed (in principle) since the fermions remain free, and one need only com- 
pute the change in the fermion two-point function. In particular, for time-independent 
perturbations we need only compute the change in the resolvent |Tf 



(7 ' 3) 

To be more explicit we must realize that the operators J X s are ill-defined due to ultraviolet 
worldsheet divergences [|14| , [T5[| so we consider matrix-model potentials SV(X) with rapid 
falloff for A — ► oo. These may be viewed as linear combinations of special state operators. 
The new S'-matrix is obtained from the old by a simple change of the bounce factor which 
is given to lowest order by (we put A = here): 

/•OO 

8Q(z) = -7T (ip-(z,X)) 2 5V(X)dX + --- (7.4) 
Jo 

The background dependence of the scattering matrix leads to background dependence 
of the location of the singularities described in section five. The variation with background 
can be computed in the A = theory using a trick. The poles zt are bound state poles, and 
indeed at zi = — i(2£ + 1) the bound state wavefunctions are just given by the continuation 
of harmonic oscillator wavefunctions: 

<MA) = H 2£+1 (2- 1 / 2 e i ^ 4 X)e- ix2 / 4 (7.5) 

which satisfy the identity: 

r°° i 

/ ^(A)?MA)dA = 2 2£+ 2 v ^(2£+l)!e-" r/4 (^ (7.6) 
Jo 
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The meaning of the ellipsis in ( [7.2[) and the hat notation is explained in [44.15]. 
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and are thus orthogonal without complex conjugation. It follows that if SV(X) has an 
expansion in terms of even powers of A (these correspond to the nonredundant operators 
2^,|TS[) then the variation of the bound-state pole with 5V is given by the expectation 



value in the normal harmonic oscillator: 

5z e = (2£ + l\6V(V2e- iw ^x)\2£ + 1) (7.7) 

where \n) is an orthonormal basis for the standard harmonic oscillator. Thus the locations 
of the singularities Z£ change as we "turn on" the special state operators displaying the 
background dependence of these singularities. 



8. 2D String Physics at High Energies 

The importance of studying string theory in the limit of ultrahigh energies has been 
emphasized by several authors [|55]-[59]]. In this section we will throw some light - literally 



- on the nature of the tachyon condensate wall by studying the reaction of the wall long 



after a single boson of energy u has been thrown at it. It was pointed out in |14j] that 
at large energies u ^> fi, with fixed string coupling, the asymptotic expansions of string 
perturbation theory cease to make sense and there must be new physics in this regime. 
We emphasize that this large order behavior of perturbation theory is in principle deriv- 
able from the continuum Liouville approach and therefore a true string-theoretic effect, 
independent of the matrix model regularization. On the other hand, the particular kind 
of new physics that arises is probably sensitive to the nonperturbative definition of the 
model. In what follows we will attempt to study this new physics concentrating on the 
nonperturbative definition we have labeled theory I with A = 0. 

8.1. Particle Production 

The unitarity equations for 1 — > 1 scattering read: 

n 

u)\S(u\ui)\ 2 (8.1) 



U z U In - LJ - Us 

n>2 JU i=l 1 i=l 



Thus, the probability that a single incoming boson of energy u produce two or more 
outgoing particles is simply expressed as 

2 



P(u) = 1 



1 

^ e *e(t«+/i)+*e((i-t)w-/i) 



o 



(8.2) 



At large u this behaves like P(u) ~ 1 — Nevertheless there is an intricate resonant 



structure in P(u) for finite u. We have plotted this function in fig. [13 
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8.2. Energy Distributions 

We can also look explicitly at the probability that two particles of energies u>i +ui 2 = u) 
are created. The result can be described in terms of a normalized energy distribution for 

D(u 2 ,uj;li) = h ^^l^l 2 (8-3) 

where 

A =T(uj 2 ; uj, -fx) - T(uj 2 ;uj, A*) 

/W2-/i /'W2+/J (8 4) 

dxe ie(u-x)+ie(x) _ / dxe ie(u>-x)+ie(x) v ' > 

-\i J fi 

Again this distribution has very interesting large-energy behavior. One finds that the phase 

integral defining T(x; u, fx) has no stationary phase point for x < \uj — a*, the integral being 

of order 1/logu. At x = hu — fx the integral turns on rapidly and is of order (ttuj/2) 1 / 2 . 

Thus we may approximate T by a step function. For — fx < u 2 < \u) + fx, A is flat and 

of order (ttuj/2) 1 / 2 . 

Outside of this range one may approximate (here we take fx <C uj 2 <C uj): 



sin(Q(fx) + Lilogui) 
A ~ — 



iu>2 logic iQ (u>2 ) 



(8.5) 



logu) 

up to some irrelevant phases. The rapidly oscillating second term leads to resonant behav- 



ior. A typical energy distribution is shown in fig. [14 

The above analysis leads to a simple approximate expression for the matrix element 
for particle creation u — > ui\ + ■ ■ ■ + uj n in the case where no u>i is small: 



S * e ™&)j™ J2 e(uj(S)-^ + f x)(-l)^ (8.6) 

ScS+ 

Note that this is much softer than the behavior obtained at any order of perturbation 
theory. 

8.3. Particle Number Distributions 

At low energies and weak string coupling it is easy to estimate the particle number 
distribution defined by the n th term in the series (|8.1| ): 

r,/ \ 2 / u\ 2n ((n - 2)!) 2 2 fuj\ 2n 1 , n . 
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This is almost a Poisson distribution: at low energies the fermions bounce with small 
phase shifts, and the subsequent bosonization proceeds via independent processes. The 
high-energy distribution function would be extremely interesting to work out. It is rather 
difficult since the approximation (|8.6Q breaks down near the boundaries of the simplex 
^ uji = u and these boundaries contribute significantly to the answer. We must leave the 
analysis of this point for the future. 

9. Conclusions 

Finally, we discuss some points of possible relevance to future work. 

The nonperturbative violation of unitarity in theory II deserves to be understood 
better. From the discussion in sec. 4.2 it is clear that theory II can be made unitary by the 
inclusion of all the soliton sectors. This suggests the following interesting open question. 
Let us begin with theory II and consider the most general way of adding degrees of freedom 
to restore unitarity. Will we be uniquely led to the free fermions of the matrix model? It 
is natural to wonder if a similar situation will some day occur in higher dimensional string 
theories, and that nonperturbative unitarity will unveil the correct degrees of freedom with 
which we should describe the theory. 

We have made some preliminary remarks on backgrounds, the main point being that 
all background dependence on V(X) can be summarized in the bounce factor e % ®( x '> v ) ^ which 
is essentially the scattering matrix of a free fermion in the potential V. It is thus clear 
that inverse scattering theory is well-suited to this problem, and we hope to investigate 
this point more thoroughly in the future.0 

Several other points must be clarified before we can achieve a complete understanding 
of background dependence. In general, the relation between the Liouville coordinate and 
the matrix model eigenvalue should be better understood. Also, the important underlying 
Woo symmetry of the c = 1 theory has played a minor role in this paper. We expect that 
a complete understanding of the background dependence will require some nontrivial use 
of this symmetry. That in turn will require a better understanding of how acts on S 
than was obtained in section seven. 

11 In this connection C. Vafa has made the interesting suggestion that the bounce factor be 
regarded as the KP tau function and that it should satisfy some difference equations analogous 
to the string equation. One could then guess that background dependence would essentially be 
KP flow on the bounce factor. 



30 



We have discovered nonperturbative dependence on parameters like the position of 
the infinite wall at A = A. Recall that in double-scaling the theory the wall is at a distance 
of order 1/N from the quadratic maximum so we are modifying the potential in the scal- 
ing region and cannot really expect physics to be independent of A. Thus, we should not 
pessimistically call the A-dependence a nonperturbative violation of universality, rather, 
we should optimistically term it the discovery of nonperturbative parameters analogous to 
the theta parameter of QCD.0 The existence of such parameters, and indeed the existence 
of different "consistent" backgrounds raises once again the persistent spectre of the unpre- 
dictability of string theory There are two possible ways this situation can be remedied. 
First, there might be further consistency conditions not yet considered which might rule 
out some possibilities. The absence-of-real-poles constraint on solutions to Painleve I is 
a case in point. In our case we have seen that the S'-matrix has a very intricate analytic 
structure. Accordingly there might be physically motivated constraints on this singularity 
structure which exclude some backgrounds from consideration. Second, some dynamical 
principle might favor some backgrounds and exclude others. Of course, this is an an- 
cient hope, but with the advent of solvable string models and a deeper understanding of 
backgrounds an answer might be forthcoming. 
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Appendix A. The Resolvent of the Upside-down Oscillator 

We list here some relevant formulae for working out the large A asymptotics of 
/(9,Ai,A 2 ) = (I(-q,Xi,X 2 ))* = i(X 1 \ w ^ = ^\X 2 ) defined in (p|). 

12 History repeats itself here: see the conclusions to the paper of Douglas and Shenker 
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This resolvent can itself be expressed in terms of parabolic cylinder functions at a 
complex energy z = [i + iq. The asymptotics of the even and odd parabolic cylinder 
functions, defined in terms of degenerate hypergeometric functions by: 



1 



^(l + e 2 ™) 1 ^ 
1 

e -i7r/8 



(W(a,x) + W(a,-x)) 



2 l/4 



r(l/4 + m/2) 



r(3/4 + ia/2) 



1/2 



e - J x 2 /4 ii?i(1/4 _ - a/2; 1/2; ^ 2 /2) 



27T 



V> (a,x) = 



e -W4| r(1/4 + za / 2 )|^M m/2) _ 1/4 (zx 2 /2) 
(W r (a,x)-W r (a,-x)) 



(A.l) 



^(l + e 2 ™) 1 ^ 
1 



^(l + e 2 ™) 1 ^ 

g — 3i7r/8 



2 3/4 



r(3/4 + ia/2) 



1/2 



r(l/4 + m/2) 

__ e -a,r/4| r( 3 /4 + - a/2) I X M . a/2A/A{ix 2 /2) 
"7T 



xe- lx2/4 iFi(3/4 - za/2; 3/2; ix 2 /2) 



(A.2) 



are most simply expressed in terms of left and right-mover "plane wave" combinations: 



— A^+oo 



T2i 



,±iF(\,z) 



_ A^ — oo 



(27rA v / l + e 2 - 2 ) 1 /2 
(27t\X\VTT^ fi ) 1 / 2> 



,±iF(\\\,z) 



(A.3) 



where 



F(X, z) = ^A 2 - zlogX + §(z) 



= \ + \ l °9 



ra+iz) 



(A.4) 



k(z) = y/l + e 2 ^ 2 - e« z 
k(z)- 1 = y/l + e 2 - 2 + . 

By demanding that the resolvent be properly normalized and vanish at infinity we see 
that if we act on L 2 (IR) then we have 



1 7T 

(Ai| — A 2 ) = -77^(Ai - A 2 )x jR (^, Ai)x L (^, A 2 ) + Ai <-> A 2 

tl — Z A 



(A.5) 
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for Im(z) > 0. From this and the above asymptotics of \ we obtain ( |2.9| ). If instead we 
choose a semi-infinite space L 2 [A, oo) we have for Im(z) > 



\h) = -0(z)9(X 1 - X 2 ) X r(z, Xi) 



if> (z,A) + 



+ Ai <-> A: 



(A.6) 



Normalizing 0(z) and using the large A asymptotics leads to the reflection coefficient: 



/r(J + 



(A.7) 



VP 



where z = fx + iq. 

In particular, for A = we recover the simpler expression 



Appendix B. An Identity on Gamma Functions 



In [14] the four point function was calculated as the limit at small loop lengths of 



the macroscopic loop amplitude. The formalism of [[14| allows calculations to all orders in 
the genus expansion. Comparing our result ( |3.13| ) to eqn. (4.38) of |b§,we find that the 
first two terms of each (corresponding to the two-bounce contribution) are equal to within 
corrections that vanish at any order. Requiring agreement of the remainder of the two 
formulas for the amplitude leads to the rather remarkable result 



n=0 



■l) n T(- qi + n) fT(q 3 + n) T(q 4 + n) 



TV. 



X 



n-qi 

ill + S 



n) 



r(?3) r(^ 4 

£(| -iLL + qi 

r(i-i/x- 




Iml e 



iirS/2 



iH + q 2 ) r(| 



ill - q 4j 



+ 



ra-ifi + q 2 ) m 



ra-iLi + q 2 + q 3 ) T(\-iLi) T(\-iLi + q 2 + q A ) T{\-in) 



(B.l) 

where S = \ Yl\qi\- This equality holds to all orders in an asymptotic expansion at large 
fj,. Expanding both sides and equating the coefficients, which are expressed in terms of 
polygamma functions, will lead to identities on these. 
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Appendix C. Small Energy and Topological Expansions 

We write here the small energy expansions of some scattering amplitudes to all orders 
in the - expansion. These are useful for some explicit unitarity checks to all orders in 
the genus expansion, since higher n-point functions enter the unitarity equations at higher 
powers of the energy. We give here the two-point function to (D(q 9 ), the three-point 
function to 0(q 6 ), the four-point function to 0(q 7 ) in two kinematic regimes and the n- 
point function to 0(q n+1 ) in the 1 — > n regime. We also give the two-point function at 
genus one, two and three. 

The coefficients at each order in the small energy expansion are written in terms of 
ip n where 

n=l ZU (C.l) 
Vn = (^)>0 



C.l. Two Point Function 

(a) Small Energy Expansion 



fJ, q R(q; -q) =qe 



8 ( 1 / 4 - 23 I 2 - 19 1 
+ Q U92(T 1 ~ 120960 ^ 3 ~ 60480 ^ 4 ~ 672(T ^ 

—29 41 1 1 

+ ^ 181440^ - 60480^^ 3 " 5040^ + 1814400 ^ + ' ' 

(C.2) 

where q is the euclidean momentum and q > 0. 
(b) Genus Expansion 

The genus one and two results are written here in a form which emphasizes that our 
results can be expanded to give integrals of n vertex operators with any choice of charges 
over moduli spaces of riemann surfaces of arbitrarily high genus. The answer is always 
some polynomial in the charges qi with coefficients that are a combination of Bernoulli 
numbers. 
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/ < T q T. q >= -^q 2 (q ~l)(q 2 -q- 1) 
J Mi, 2 24 



(C.3) 



f 2 ^ 

/ < T q T_ q >= -|— TT (g - r) (3g 4 - lOg 3 - 5q 2 + 12q + 7) 

JM 2 ,2 5760 r =l 



(C.4) 



M 3 , 



< 99 > 2903040 



Y[(q - r)(9q 6 - 63q 5 + 42q 4 + 217q s - 205q - 93) 



r=l 



(C.5) 



C.2. Three Point Function 



fi^R( qi ,q 2 ;q 3 ) 



e lqs ^°qiq 2 \q 3 \ 



^i - y^(9i + 9i92 + gl)^3- 



12 



|g 3 |(2^ + <?i<? 2 + q\)^\^2 + 



(C.6) 



where 9i, 92 > and 53 < 0. 



C.5. Four Point Function 



The four point function in the 1 — > 3 and 2 — > 2 kinematic regimes is given up to 
0(q 7 ) where we can check that one particle irreducible amplitudes of [|39| are not analytic 
at genus one. 



fi M R( qi , 92, q 3 ; 94) =e^ qi q 2 qM 



^2 + \qM - 57(9? + q\ + qj + 94)^4- 



2 1 3 2 2 3 2 2 

9293 9i9| 9293 , 9jh_ 

2 2 2 3 J 
^1^3 , 9i 5gfg 2 5gig% g| 5gfg 2 5g^g 3 

2 1 2 6 6 2 6 6 

„ „ „ 1 5 gjg3 5gigj 5g 2 gj gg, 

919293 H ^ 1 7 1 7. r —) + • • 

O O O 2 



(C.7) 



In the following q\ = max \qi\ 
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^+i>R(q 1 ,q 2] q 3 ,q i ) 



--e^° qi q 2 \q 3 \\q 4 \ 



- qi^r 



12 



4/„2 



(Qi + Q1Q2 + q 2 + Q1Q3 + 9293 + q 3 , 



,2/9i , Q1Q2 , 9i 92 , 92 , 9?93 , 919293 9i93 
2 6 4 4 12 4 4 4 

4^ 4 ^ 4 T i2 4 J 



(C.8) 



C.^. N-point Function 
For n > 2 we have: 



/i l9n| i?(gi, g 2 , • • • ! 9n) =9l92 • • • gn-ik- 

1 1 n-2 



n-2 



2 •-«— ' \ r 

r=0 



E ( n - 2 ) V'n-2-rV'r) + ' " 



(C.9) 



These first two terms can easily be obtained by using the recursion relation 



R(qi,q 2 ,...;q n ) ~ qi —R(q 2j . . .■ q n ) 
for — > 0, and they are precisely what we expect from unitarity. 



(CIO) 



Appendix D. Proof of Equivalence to Bosonization 

We wish to prove that the S'-matrix we compute may be obtained from the simple 
scattering amplitude ( |4.4| ) for free fermions via the bosonization prescription of Q4.2| ). We 
will construct the proof in Euclidean space to connect to the derivation of the 'filtration for- 
mula'. As with that formula, all results are later continued to Minkowski space. Formally, 
we need to show 



n d ^ n a ^ + ^) at + & + 9i)ii"> 



ti k n k 

E E Ti/ ^ j f-(T j ,-( j )H^^ + C j + q(T j ))a^ + ( j )\ f ,) 
k=iT 1 n...uT k =s- j=x 3=1 



(D.l) 
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where S = {q±, . . .q n }. Note that since f-(T,Q) oc —9(—Q)9(Q — the left-hand 

side of ( p.l| ) is normal-ordered. The proof proceeds by induction on n. For n = 1, the 



claim is explicitly verified in (|4.2|) . Assuming the assertion holds for n — 1 tachyons, we 
proceed to n by multiplying the "fermionized" state corresponding to S = {q2, ■ ■ - qn} on 
the left by J d£xa(fi + £i)a'(n + £1 + q\). To do this, we split the integral into a normal- 
ordered part 1 1 = f Q qi d£i and the remainder I2. In I±, the operators can be commuted 
to the right to canonically order the result, leading to 

71— 1 » k 

h=Y, e / n^n^^-^H^+^+^w^^) i/*> ■ ( D - 2 ) 

k=i Tin...nr fc =S" i j =1 
Tfe={gi} 

In I2, the operators on the left annihilate so the contribution is a sum of anticommu- 
tator terms. Using the identity 



/_ (T, Q) = /_ (T \ {g* } , Q - q* ) - f_ (T \ {q* } , Q) 



(D.3) 



for q* G T, one verifies that the jth contribution is 

i>oo \ 

+ / a{fi + ^ 1 )a^(fi + ^ 1 + q 1 ),a^(fi + $ j +q(T j ))a(fi + ^) 

-00 J—qi/ 

Adding the two contributions proves the claim by induction. 



;d.4) 
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Figure Captions 

Fig. 1. a.) A pictorial version of the integral / for positive momentum, (b) A pictorial 
version of the integral / for negative momentum 

Fig. 2. Incoming and outgoing vertices. The dotted line carrying negative (positive) 
momentum qi should be thought of as an incoming (outgoing) boson with energy 
\qi\. Momentum carried by lines is always conserved as time flows upwards. 

Fig. 3. The possible real histories for 1 — > 2 scattering 

Fig. 4. Some possible real histories for 2 — > n scattering 

Fig. 5. Diagrammatic rules in Euclidean space 

Fig. 6. Diagrammatic rules in Minkowski space 

Fig. 7. Typical bracket configurations 

Fig. 8. A real history as a composition of three maps 

Fig. 9. An incoming boson dissociates 

Fig. 10. Analytic structure in the complex s plane 

Fig. 11. A set of branch cuts for S(u)\lj), illustrated here fore the case A = 0. 

Fig. 12. Two paths defining the 1 — > 1 amplitude on the second and third sheets. 

Fig. 13. The amplitude 1 — P{ui) for fi = 8 as a function of u>. The complicated structure 
observed here is probably not an artifact of the existence of the wall. Indeed we 
may also plot the amplitudes in theory II to obtain similar pictures. 

Fig. 14. Energy distribution integral \ A\ 2 for fi = 5 and u> = 6. Again, similar results hold 
in theory II. 
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